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We investigate the behavior of the chiral condensate in lattice QCD at finite temperature and 
finite chemical potential. The study was done using two flavors of light quarks and with a series 
of [3 and ma at the lattice size 24 x 12^ x 6. The calculation was done in the Taylar expansion 
formalism. We are able to calculate the first and second order derivatives of ('ip'ip) in both isoscalar 
and isovector channels. With the first derivatives being small, we find that the second derivatives 
are sizable close to the phase transition and the magnitude of ip'ip decreases under the influence of 
finite chemical potential in both channels. 


I. INTRODUCTION 


Since the advent of lattice QCD in the seventies by 
K.G. Wilson [1], it has been proved that the theory is ex¬ 
tremely successful in the analysis of low-energy dynam¬ 
ics among mesons and baryons [2 E]. Combined with 
large-scaled computations on supercomputers, people 
have been investigating various non-perturbative quali¬ 
ties, such as hadron spectrum, chiral transitions, behav¬ 
ior of the glueballs, hadronic matrix elements, spatial mo¬ 
mentum dependence and vector current correlation[4l|5]. 

However, lattice QCD also suffers from some shortcom¬ 
ings. Eor example, it violates some of the important sym¬ 
metries that the continuum theory acquires and which 
can only be restored in the continuum limit. One of the 
most important symmetry is the chiral symmetry. A well 
known no-go theorem due to Nielson and Ninomiya [6] 
says that chiral symmetry has to be realized differently 
on the lattice. The most elegant way is known and the 
fermion matrix of the lattice theory has to satisfy the 
so-called Ginsparg-Wilson relation [7]. One example in 
this category is the so-called overlap fermion. However, 
practical simulations of overlap fermion encounters other 
technical problems and is rather costly. A compromise to 
this problem is to use staggered fermion which preserves 
part of the continuum chiral symmetry [8] and runs ef¬ 
fectively on supercomputers. 

In this letter, using staggered quarks, we study the 
chiral condensate {'ipip) and the relevant chiral symmetry 
breaking. As is well-known, is an important order 

parameter in the measurement of phase transition in lat¬ 
tice QCD [9llTT|. Basically, the quantity {'ip'ip) exhibits 
a fast decrease around the critical point but placidity at 
other place. This enables us to investigate the phase tran¬ 
sition and chiral properties in both low temperature and 


QGP phase. Simulations of lattice QCD at finite density 
encounters another well-known long-lasting problem-the 
sign problem-the name borrowed from condensed matter 
physics where it also appears in the simulation of models 
with fermions. In this letter we will follow the strategy 
of Taylor expansion method [12] in which all physical 
quantities are expanded around /i = 0. 

This paper is organized as follows. In the next section, 
we briefly review the derivation of the derivatives for the 
chiral condensate. In the third part, our numerical re¬ 
sults for the chiral condensate are presented and we will 
summarize and conclude in the fourth part. 


II. TAYLOR EXPANSION 


When the chemical potential is present, the simula¬ 
tion encounters the infamous sign problem. However, if 
the chemical potential is not too large, everything can 
be expanded into a Taylor expansion in (yu/T). Eor the 
purpose of this letter, we need the expansion for ('^'0), 
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The first order and second order responses in the above 
expansion can be computed using numerical simulation 
at vanishing chemical potential. Using this equation, we 
can investigate the behavior of {'ipip) at a small but non¬ 
vanishing chemical potential. 


1 












2 


The staggered quark fermion matrix is given by, 

D{U, fi) = maSn,m + ^ Vain) 

a=x,y,z 

^ ^)^n—a,m] 

+ \vt{n)[Ui{n)e^5^^f „^ - Uj{n - i) 

where m is the bare quark mass and 77 ^^^ is a parameter 
only depend on the parity. The parameter jl = ji/^NtT) 
designates the chemical potential. It is switched on for 
each flavor of quarks. In our simulation, we have both u 
and d quarks whose chemical potential will be denoted by 
jj^u and iJLd respectively. The quark propagator is related 
to the Dirac operator D\U]jl] in the background gauge 
field configuration U as: 

g{ii) = D{U,ii)-^ (3) 

In lattice QCD, chiral condensate can be written as fol¬ 
lows: 

= Re{G) = (^^Re{TT[g{M] + lV[5(Ad)]}^ , (4) 


where Tr implies summing over all indices and (•) indi¬ 
cates averaging over the gauge field ensembles. 

For the observable (G), one can take the first and 
the second order derivatives with respect to the reduced 
chemical potential jl = ji/T and obtain 
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where A is the fermion determinant given by 
A = det(D(/i„)) det(D(/id)) , 


( 5 ) 


(6) 


and a dot indicates the derivative with respect to jl. 

It is interesting to investigate the response of the chiral 
condensate to both the isoscalar chemical potential (the 
traditional finite density case) and its isovector counter¬ 
part (the finite isospin density case). For the isoscalar 
we set 


jds = idu = jdd , (7) 

while for the isovector case we choose 


jdv = idu = —jdd • 


( 8 ) 


We will call these two cases as isoscalar and isovector 
channel in the following by calculating the corresponding 
derivatives of G with respect to the corresponding chem¬ 
ical potential. It is known that the first order derivative 
for both the isoscalar case and the isovector case are zero, 
so the second order derivatives become crucial in our cal¬ 
culation. 
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TABLE I. Simulation parameters used in this study. All lat¬ 
tices are of the size 24^ x 12 x 6 and Nconf stands for the 
number of configurations. 


ma /? 

0.020 5.292 0.367(13) -1.71(12) -1.96(18) 

0.020 5.327 0.351(13) -2.26(31) -2.75(30) 

0.020 5.347 0.325(28) -5.73(42) -6.21(34) 

0.020 5.373 0.138(24) -6.73(45) -7.12(53) 

0.020 5.392 0.129(10) -2.27(19) -2.25(22) 

0.020 5.422 0.119(8) -1.43(5) -0.98(4) 

0.015 5.317 0.367(26) -2.81(31) -2.79(29) 

0.015 5.337 0.335(33) -7.07(59) -6.55(56) 

0.015 5.363 0.144(28) -6.39(36) -5.92(34) 

0.015 5.382 0.123(8) -2.13(20) -2.62(24) 

0.012 5.327 0.377(21) -2.97(26) -2.83(28) 

0.012 5.347 0.347(34) -6.08(41) -6.17(45) 

0.012 5.373 0.148(28) -5.68(35) -5.67(32) 

0.012 5.392 0.126(7) -2.55(23) -1.84(17) 


TABLE 11. The values of {'ij'ij) and its second order derivatives 
in the isoscalar and isovector channel. 


III. NUMERICAL SIMULATION AND RESULTS 


The conventional Wilson plaquette action is used 
which is characterized by the parameter The gauge 
field configurations are generated using conventional R- 
algorithm of Hybrid Monte Carlo with molecular dynam¬ 
ics step ^ = 0.01 and trajectory length of 50. The size of 
the lattices are all 24 x 12^ x 6 and for each parameter set, 
600 gauge field configurations are obtained. By scanning 
the temperature dependence of the Polyakov loop the 
ratio T/Tc could be determined. This information, to¬ 
gether with other simulation parameters are summarized 
in Table 1. 
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We have calculated the value of and its second 

order derivatives to both fls and fiy foi* all our data sets 
and the results are listed in Table 11. Temperature depen¬ 
dence of the second order derivatives of {'ip'ip) with respect 
to fls and fly are illustrated in Figure 1 at ma = 0.020. 
It is seen from the figure that the second order deriva- 



FIG. 1. Second order derivatives of at ma = 0.020. 


tives in both isoscalar and isovector channel share some 
common features. They are both negative but the abso¬ 
lute value is very small when away from the critical point. 
Close to the critical point, both increase substantially. 

Now it is possible to consider the behavior of {ipip) 
with chemical potential and temperature in the critical 
region for both channels. We omit the first order deriva¬ 
tive for the quantity since they vanish and only consider 
the second order derivatives. Using the data listed in Ta¬ 
ble II, we obtain the following expansionw near pc- For 
instance, slightly below Tc oX P = 5.347 and ma = 0.020, 
we obtain: 
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As we see, the derivative corrections are quite substantial 
in the critical region. Thus, the effect of the chemical 


potential makes Pc to drop from its original value at /r = 
0. At an even lower temperature, for example, p = 5.292 
and ma = 0.020, the data suggests 


^ 79 2(2.7) - 5.13(37) 

(thyiiiv) 
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These derivative corrections are not that large as com¬ 
pared to the case in the critical region. In the phase above 
Tc (in the QGP phase) at P = 5.422 and ma = 0.020, we 
obtain 
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25.8(1.8) - 2.95(12) 
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We can now plot the results for the chiral condensate 
at small chemical potential in both the isoscalar and the 
isovector channel. This is illustrated in Figure 2. In this 
figure, we include the behavior of for the chemical 
potential ps^v = 0.05 and ps,v = 0.10. Since the sec¬ 
ond order derivatives are all negative, we find the crit¬ 
ical temperature tends to decrease under the influence 
of ps,v- Thus, in the low temperature phase, turning 
on the chemical potential brings the system closer to the 
phase transition where chiral symmetry is restored and 
decreases the magnitude of the chiral condensate. In the 
high temperature phase, however, since chiral symmetry 
is already restored, the responses of the chiral conden¬ 
sate to the isoscalar and isovector chemical potential are 
relatively small. 


IV. CONCLUSIONS 

In this work, we have studied the response of the chi¬ 
ral condensate ('^'0) to the chemical potentials using the 
Taylor expansion method. The quantity is expanded 
around /x = 0 and the second derivatives of {'ip'ip) with re¬ 
spect to both fls and fiy are obtained. As is seen, though 
the first order derivatives are small, the second order re¬ 
sponses are sizable and exhibit several unique features. 
The behavior of the responses for the ('^'0) is closely re¬ 
lated to the chiral restoration. For both isoscalar and 
isovector channel, we find that the critical temperature 
decreases with the influence of both chemical potentials, 
and the magnitude of ('^'0) tends to decrease under the 
finite chemical potential. 
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FIG. 2. Behavior of at finite isoscalar chemical potential 
(top panel) and the isovector potential (bottom panel). 
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